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PREFACE 

Part  of  the  research  program  of  The  RAND  Corporation 
consists  of  basic  supporting  studies  in  mathematics.  This 
Memorandum  presents  a  method  for  handling  a  problem  that 
frequently  arises  in  mathematical  physics  in  connection 
with  equations  of  heat,  radiative  transfer,  etc.,  namely, 
that  of  replacing  a  system  of  nonlinear  functional 
equations  with  one  that  is  linear- 


One  of  the  basic  problems  of  mathematical  physics  is 
that  of  replacing  a  nonlinear  functional  equation  by  a 
more  tractable  (analytically  and  computationally''  linear 
equation.  More  generally,  one  wants  to  replace  a  system 
of  nonlinear  equations,  often  derived  from  a  single 
equation  by  means  of  an  expansion  in  a  parameter,  an 
orthogonal  expans ionT-^th  a  system  of  linear  functional 
e<]uations.  To  treat  this  ^losure^  problem,  the  authors 


present  a  method  basea  upon  the  concept  of  relative 
invariants  and  the  use  of  the  multidimensional  Lagrange 


expansion  theorem. 
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RELATIVE  INVARIANTS  AND  CLOSURE 


I,  INTRODUCTION 

One  of  the  basic  problems  of  mathematical  physics  is 
that  of  replacing  a  nonlinear  functional  equation  by  an 
analytically  and  computationally  more  tractable  linear 
equation.  More  generally,  the  problem  is  that  of 
replacing  a  system  of  nonlinear  functional  equations  by  a 
system  of  linear  functional  equations.  Often  this  system 
of  nonlinear  equations  if;  derived  from  a  single  equation 
by  means  of  an  expansion  in  a  parameter,  or  of  an  ortho¬ 
gonal  expansion. 

Consider,  for  example,  the  nonlinear  heat  equation 

with  the  boundary  and  initial  conditions 
(1.2)  u(x,0)  “  g(x),  0  <  X  <  1, 

u(0,t)  »  u(l,t)  “0,  t  >  0. 


If  we  set 

oo 

(1.3)  u  =  E  u  (t)  sin  ntrx, 

1  n 
n*l 


substitution  in  (I.l)  yields  the  infinite  system  of 
nonlinear  equations 


u;(t) 


—  n  u  + 
n 


00 

E 
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a,,  u,u. 
ijn  i  j 


u  (0) 
n'  ' 


n 


(1.4) 


•  •  • 
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^2 


where 

CD 

g(x)  'v.  r  g  sin  nTTx, 
n-1  ^ 


We  now  faee  the  problem  of  closure  (see  [1],  [2]), 
which  is  to  say  of  replacing  (1.4)  by  a  finite  system  of 
ecjuations,  preferably  linear.  We  shall  present  a  method 
based  upon  the  concept  of  relative  invariants  and  the  use 
of  the  multidimensional  Lagrange  expansion  theorem  [2], 

[3].  In  some  cases,  as  in  the  case  abcjve,  a  preliminary 
use  of  approximation  methods  may  be  required, 

2.  FINITB-DIMENSIONAL  CASE 

In  order  to  illustrate  the  idea  of  relative 
invariants,  let  us  first  consider  a  finite  system  of 
nonlinear  differential  ecju^^tions,  say 

(2.1)  ^  +  g^(x),  Xj^(O)  -  c^,  i  =  1,2,...,N, 

where  we  suppose  that 

(2.2)  (a)  Re(xp  >  0, 

(b)  |c^|  is  sufficiently  small, 

(c)  gj^(x)  is  a  power  series  in  the  components  of 
X  beginning  with  cjuadratic  terms. 

We  shall  impose  another  condition  upon  the  X^  in  a 


moment. 
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Our  aim  is  to  find  N  functions  of  the 
{(p^(x)l,  with  the  properties  that 

dq). 

(2.3)  ST  --Vi- 
and 

(2.4)  q)^(x)  =  x^  +  h^(x), 

where  the  h^(x)  are  power  series  in  the  components  of 
X  beginning  with  quadratic  terms. 

If  such  functions  exist,  we  have,  from  (2.3), 

-X.t 

(2.5)  cpj^(x)  “  e  qp£(c),  1.  ~  X,2,*..,N. 

Hence,  the  x^  are  obtained  as  solutions  of  these  N 
simultaneous  equations.  Since  from  (2.4)  these  equations 
have  the  form 

-X.t 

(2.6)  =  e  ^  (p^(c)  -  h^(x),  i  =  1,2,.. .,N, 

we  are  in  an  ideal  position  to  employ  the  raultidiraensiontil 

Lagrange  expansion  to  obtain  a  power  series  expansion  for 

-X.t 

the  in  terms  of  the  quantities  e  cp^(c.).  Cutting 

this  off  at  any  desired  stage,  we  have  the  required 
closure. 

Furthermore,  if  we  so  wish,  we  can  employ  the 
Lagrange  expansion  theorem  to  calculate  N  additional 
functions  ^^(x),  i  “  1,2,...,N,  in  place  of  the  N 

^N' 


functions  Xj^,X2i.../ 


3.  WHAT  1$  REQUIRED  FOR  THIS  FOKMA1.ISM? 

We  mentioied  above  that  a  further  condition  upon  the 
was  required.  In  order  to  c.ee  how  this  arises,  let 
us  examine  the  problem  of  determining  the  functions  q)j_(x) 
of  (2.3)..  Write 

2  2 

(3.1)  cp^(x)  "  +  2a^2i^l^2  ^221^2  +  •••, 

and  let  us  restrict  ourselves  to  the  case  i  =  1.  We 
have 

(3.2)  "at"''  ar”  ^®i2i  3t“  ^2 

dx2  dx- 

^^121^1  W~  ^^221^2  IT* 

“  "  H^l  ^111^1  ^12i’^l’^2 

+  X^x^  +  giiiy^i  +  ^121^1^2  * 

+  2a|^2]^^2^~  ^1^1  ^111^1  ^121^1^2 

+  2aj^2i^lt~  ^2’'2  ®112’'l  ^ 

+  2a22]^5C2[—  X2X2  +  •**]• 

L-^t  us  concentrate  solely  upon  first  and  second  order 
terms,  and,  in  partici’lar,  upon  the  determination  of 
^121^  and  ^221*  have,  from  (3.2), 

<Pl(x)  “  -  +  g]^2l’'l’'2  8221’‘2  ■*■ 


(3.3) 


-5- 


2  2 

There  are  do  other  terms  involving  ^2' 

If  (3.3)  is  to  hold,  for  i  “  1  we  imist  ha/e,  upon 
equating  coefficients, 

(3.4)  =  -  ^1^111^ 

gl2i  ”  ('-^A^  +  2X2)aj^2i  “  H^121" 

§221  “  ^^2^221  “  “  ^1^221* 

These  relations  determine  the  quantities  ^221 

uniquely,  provided  that 

(3.5)  jt  0,  X^  +  2X2  ^  0,  2X2  -  X^  0. 

As  far  as  the  determination  of  all  of  the  coefficients  is 
concerned,  we  see  that  we  have  the  condition 


(3.6) 


'^l^l  ^2^2 


+  0' 


for  any  nontrivial  set  of  positive  or  negative  integral 
or  zero  values  of  the  m^.  This  is  a  frequently  met 
condition. 

This  condition  is  not  merely  a  reflection  of  the 
weakness  of  the  method.  If,  for  example. 


(3.7)  ^  ~  ^1^  x^(0)  =  Cj^, 

X2  =  —  2X2  +  Xp  X2(0) 

-'2t 

we  see  that  a  term  of  the  tc 
x^  and  X2 


i. 

enters. 

t 


The  functions 


and 


are  not  power  series  in  e 


Observe  that  we  face  this  situation  in  treating  (1.4), 
2 

since  here  X  *  n  X, . 

n  1 

4 .  PRELIMINARY  CLOSURE 

To  avoid  a  violation  of  the  condition  in  (3.6),  we 

can  use  some  preliminary  closure  techniques.  Suppose, 

for  example,  that  2X2  —  X^  =  0,  which  means  that  the 

third  equation  in  (3.4)  is  troublesome  unless  S22I 

^221  both  zero.  If  §221  ~  very  simply 

set  ^221  *  ^221  ^  approximation 

2 

technique  to  replace  the  §221^2  ^  linear  combination 

of  other  terms  appearing  in  the  expansion  of  g|^(x); 
see  [1], 

Thus,  we  can  set 
(^•1-)  ®221^2  *  ^11^1  ^12^2 

where  the  c. .  are  to  be  determined  by  a  least— squares 
fit.  Various  self-consistent  methods  can  now  be  used. 

5.  INFINITE-DIMENSIONAL  CASE 

Consider  now  the  infinite  system  of  ordinairy 
differential  equations 
dx . 

(5.1)  ^  =  -  X^x^  +  g^(x),  x^(u)  »  c^,  i  *=  1,2, 

obtained  from  a  nonlinear  operator  equation 
Uj,  ■=  Au  +  g(u). 


(5.2) 


-1 


where  we  have  used  an  expansion  for  u  which  diagonalizes 
the  operator  A. 

It  is  easy  to  see  that  if  the  condition  of  (3.6)  is 
met,  we  can  determine  the  coefficients  of  the  function? 
cp^(x)  recurrently,  solving  a  single  equation  for  a  single 
unknown  at  each  step.  The  determination  of  qp^Cx)  can  be 
carried  as  far  as  desired.  The  stopping  '*ule  is  now  a 
closure  technique.  Having  determined  the  N  functions 
q)^(x),  i  -  1,2,...,N,  to  some  degree  of  accuracy,  the 
N  equations 

— X.  t 

(5.3)  qp^(x)  “  e  qp^(^)^  ^  **  1>3,...,N, 

may  be  solved  as  indicated  above,  or  in  some  other  way,  to 
yield  the  x^,  i  =*  1,2,...,N.  We  suppose  in  (5.3)  that 
on  the  left-hand  side,  all  terms  involving  x^, 
i  >  N  +  1,  have  been  omitted. 

There  are  many  variants  of  this  procedure,  each 
equivalent  to  a  different  closure  method. 

6.  APPLICATION  TO  THE  NONLINEAR  HEAT  EQUATION 

In  applying  this  method  to  (1.4),  we  see  that  the 
first  violation  of  (3.6)  occurs  in  connection  with  the 
determination  of  fourth-order  terms  in  the  functions 
qp^(x)  (since  Hence,  if  the  deviation  from 

linearity  is  small,  or  equivalently  if  max  jg(x)j  is 

X 

small,  the  method  sketched  above  will  yield  the  second— 
and  third— order  correction  terras  for  the  solution  of 


(6.1) 


sin  nTTX 


_g_ 


u,  (x,t)  =  T.  g  e 

1  rx^i  n 


in  a  convenient  fashion. 
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